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ABSTRACT 


The  purpose  of  this  study  was  to  determine  the  feasibility  of 
computing  residual  displacements  in  runways  and  crater  repairs  under 
repeated  aircraft  loadings.  This  report  presents  a  methodology  for 
computing  residual  displacements  that  is  economical  and  easy  to  use.  It 
also  contains  a  literature  review  and  a  presentation  of  a  second 
procedure  to  compute  residual  displacements.  Sample  calculations  were 
done  to  illustrate  the  first  procedure  and  the  results  compared 
favorably  with  previously  published  test  data.  The  limitations  of  the 
second  procedure  were  also  discussed. 
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CHAPTER  1 

Scope  and  Concept  for  Computing  Residual  Displacements 

Introduction 

The  purpose  of  this  report  is  to  present  the  results  of  a  feasi¬ 
bility  study  for  computing  the  residual  displacements  in  airport  pave¬ 
ments  and  crater  repairs.  The  study  emphasized  the  application  of  the 
finite  element  method  of  analysis  with  the  incorporation  of  non-linear 
behavioral  laws  to  represent  the  base  and  subbase  material  of  the 
pavement  structure.  The  report  consists  of  three  sections.  The  first 
section  presents  the  scope  and  concept  of  the  study.  Section  two 
presents  a  brief  literature  review  and  background  information.  Section 
three  is  a  presentation  of  several  proposed  procedures  along  with  sample 
calculations.  Limitations  and  further  research  for  each  procedure  are 
also  presented.  In  general,  the  report  concludes  that  there  does  exist 
an  ec  lomical  and  feasible  procedure  for  computing  residual  displace- 

N 

ments. 

Scope  and  Concept  of  a  Method  for  Computing  Residual  Displacements 
Concept 

Airport  pavements  are  subjected  to  many  different  types  of  load 
applications.  Runways  are  subjected  to  high-speed  takeoff  and  landings, 
l.e.,  dynamic  and  impact  loading  from  many  different  types  of  aircraft. 
Taxiways  handle  aircraft  under  static  loading  conditions.  High-speed 
runway  exits  handle  dynamically  applied  wheel  loadings. 

An  idealized  analysis  procedure  that  used  the  finite  element  method 


is  presented  below.  This  procedure  presupposes  that  capital  costs  do 


not  hinder  testing  of  pavement  materials  or  computer  operating  time,  and 
the  computer  core  size  does  not  limit  the  size  of  the  finite  element 
mesh.  First,  the  dynamic  and  static  constitutive  (stress-strain) 
properties  of  the  pavement  materials  would  be  determined.  Secondly,  a 
three  dimensional  finite  element  grid  could  be  constructed.  Lastly,  the 
actual  plane  movements,  including  takeoffs,  landings  and  taxiing,  could 
be  simulated  in  the  time  domain,  taking  into  account  the  fact  that 
planes  traffic  at  different  positions  along  the  length  of  the  runways 
and  along  different  positions  with  respect  to  the  center  line  of  the 
runway.  The  computer  would  actually  simulate  the  landing,  takeoffs  and 
taxiing  in  three  dimensions  with  aircraft  moving  simultaneously  through¬ 
out  various  portions  of  the  airport.  The  computer  would  automatically 
switch  from  dynamic  constitutive  relations  to  static  constitutive 
relations  depending  on  the  velocity  of  the  aircraft. 

The  hypothetical  description  of  a  mythical  computer  program  could 
be  achieved  on  perhaps  a  few  computers  in  the  country.  Such  a  program 
would  not,  however,  be  possible  to  operate  on  a  majority  of  the  main¬ 
frame  computers  at  most  government  installation  or  universities.  Such  a 
program  may  be  the  state  of  the  art  in  the  near  future,  but  this  is  not 
consistent  with  the  current  state  of  the  art  exposed  during  the  present 
research  project.  The  procedure  developed  during  the  study  was  well 
within  the  financial  limitations  and  compc.er  size  limitations  as  they 
presently  exist. 

Current  finite  element  analyses  of  airport  pavements  include 
axisymmetric ,  plane  strain,  and  prismatic  solid  representation.  Elastic 
layer  theories,  mats  on  a  Winkler  subgrade  and  visco-elastic  solutions 
are  also  currently  used.  Most  computation  procedures  handle  a  one-time 


loading  of  Che  pavement  system.  Many  procedures  use  a  resilient  modulus 
so  that  the  response  during  one  of  the  many  load  applications  can  be 
analyzed.  None  of  the  current  programs  had  the  capacity  to  handle 
multiple  application  of  cyclic  loadings. 

Pavement  materials  include  concrete,  bituminous  asphaltic  concrete, 
gravel  subgrades  and  soil  subbases.  These  materials  actually  have  a 
nonlinear  relation  between  stress  and  strain;  that  is,  when  a  load  is 
applied  and  then  removed,  the  material  will  not  return  to  its  initial 
position.  The  material  accumulates  a  permanent  strain  or  set.  A 
hypothetical  stress-strain  diagram  is  shown  in  Figure  1.  The  loading 
and  unloading  path  are  indicated  as  well  as  the  permanent  strain. 

When  a  load  is  applied  to  a  pavement  material  and  removed  many 
times,  as  occurs  due  to  traffic  applications,  the  permanent  strain 
accumulates  leaving  a  permanent  or  residual  displacement  in  the  pave¬ 
ment.  Such  ruts  are  often  visible  in  unpaved  roads  or  in  bituminous 
highway  pavements.  The  permanent  strains  usually  accumulate  rapidly 
during  the  first  few  applications  of  the  load.  After  repeated  appli¬ 
cations,  the  accumulated  set  strain  will  either  reach  a  limiting  value 
or  will  undergo  a  rapid  Increase  due  to  a  fatigue  type  failure.  An 
intermediate  case  where  the  strains  increase  according  to  a  prescribed 
function  is  also  possible.  Figure  2  illustrates  these  behavioral 
patterns.  Fatigue  failures  have  been  discovered  to  be  related  to  the 
magnitude  of  stress  application  as  compared  to  ultimate  strength  of  the 
material.  Some  highway  materials  also  exhibit  a  creep  behavior. 

In  laboratory  testing  studies,  many  investigators  have  been  able  to 
relate  the  accumulated  set  strain  to  the  set  strain  produced  during  the 
first  loading  cycle.  In  general  the  accumulated  set  strain,  er,  is  a 
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function  of  the  number  of  loading  N,  the  ratio  of  the  stress  appli¬ 
cations,  Aa,  to  the  ultimate  strength  o^ ,  confining  pressure,  o^,  and 
the  first  set  strain,  e^, 

er  -  F(e1,Ao/a1,  o3,  N)  (1) 

Compositional  changes  also  affect  the  accumulated  strain  relationships. 

The  existing  finite  element  programs  could  be  slightly  modified  to 
calculate  the  permanent  strain  after  the  first  cycle.  The  accumulated 
set  strains  for  each  element  could  then  be  calculated  using  Equation  1 
and  the  total  surface  settlement  Ah  could  be  calculated  by  summing  the 
change  of  height  of  each  element,  or 

No.  of  element 

Ah  =  Z  e  *h  (2) 

i=l 

where  e  is  the  accumulated  set  strain  for  each  element  and  h  is  the 

fj  1  1 

height  of  each  element.  It  is  proposed  that  programs  use  non-linear 
constitutive  laws  to  represent  the  material  behavior.  Linear  consti¬ 
tutive  laws  are  valid  if  appropriate  engineering  judgment  is  exercised 
in  selecting  the  secant  or  tangential  material  modulus  to  approximate 
the  non-linear  behavior.  As  will  be  shown  later,  the  initial  strain 
finite  element  procedure  can  also  be  used  to  calculate  permanent  dis¬ 
placements. 

Scope 

Equation  2  represents  the  current  concept  described  in  the  litera¬ 
ture  and  was  found  to  be  the  most  feasible  procedure  most  easily  adopt- 
able  for  computing  residual  displacements.  It  was  the  scope  of  this 
project  to  compare  the  application  of  this  concept  with  published  data 


,V’WT  • 
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regarding  rue  accumulations  and  permanent  deflections  in  bomb  crater 
repairs  and  airport  pavements.  The  use  and  the  problems  of  using  the 
initial  strain  procedure  were  also  studied. 


CHAPTER  II 


Background  Information 


Introduction 

This  section  of  the  report  provides  a  summary  of  some  of  the 
pertinent  literature  that  was  reviewed.  The  section  is  written  in  the 
chronological  order  of  the  development  of  the  various  ideas  regarding 
the  analysis  of  pavement  systems  and  the  accumulation  of  permanent 
deflections. 

Literature  Review 

Ahmed  and  Larew  (1962)  publish  the  results  of  static  and  repeated 

load  triaxial  testing  on  three  types  of  subgrade  soils.  These  soils 

were  a  limestone  residual  silty-clay,  a  sandy  clay,  and  a  micaceous 

silt.  All  their  testing  was  done  with  a  single  confining  pressure  of 

10  psi.  Duplicate  samples  were  used  for  each  soil.  The  average  peak 

devia'.or  stresses  at  static  failure  were  determined  for  each  soil,  Ao  . 

s 

A  series  of  repeated  load  tests  were  run  on  each  soil  where  the  repeated 

deviator  stress,  Ao^.,  was  held  at  certain  ratios  of  the  static  ultimate 

deviator  stress.  Studies  were  also  done  to  show  the  effects  of  water 

content,  density  and  compactive  effort. 

The  results  of  the  study  showed  that,  for  the  ratio  of  the  dynamic 

deviator  stress  to  the  ultimate  static  dewater  stress,  Aa  /A a  ,  below  a 

r  s 

critical  level,  the  deformations  reached  a  limiting  value  after  a 
certain  number  of  load  applications.  For  a  Aa^/Aag  ratio  above  a 
critical  value,  a  greater  number  of  repetitions  would  increase  the  rate 
of  deformation  per  load  application  and  the  sample  would  fail  after  a 


i 
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certain  number  of  load  applications.  Figure  3  illustrates  these 
results. 

Hardin  (1971)  reports  the  results  of  using  specially  designed 

equipment  for  determining  shear  stress  relations  of  soil.  The  1971 

report  presented  a  model  for  a  silty  sand  and  a  silty  clay.  In  1972  and 

1973,  the  model  was  applied  to  additional  soils  and  gravel. 

The  model  relates  the  sharing  stress,  t,  to  the  shearing  strain,  y, 

using  a  secant  shearing  modulus  Gg.  See  Figure  4.  Gg  was  expressed  as 

a  hyperbolic  equation  as  G  -  G  /(1+y,  )  where  G  was  the  maximum  or 
r  n  s  max  h  max 

initial  tangent  shearing  modulus,  and  <5^  was  the  hyperbolic  shear 
strain, 


_3 L  n  + _ * _ i 

Yr  L  exp(y/yr)J 


(3) 


where  6^  is  the  reference  shear  strain  and  a  is  an  experimentally 
determined  parameter. 

The  reference  shear  strain,  yf,  was  defined  as  the  ratio  of  the 

maximum  shearing  stress,  t  to  G  The  advantage  of  Hardin's  model 

max  max 

was  that  the  results  of  many  different  strength  tests  could  be  non- 
dimensionalized  in  terms  of  the  reference  strain  as  shown  in  Figure  5. 
The  model  was  related  experimentally  to  such  characteristics  as  soil 
type,  void  ratio,  plastic  index,  degree  of  saturation,  rate  of  loading, 
and  number  of  loading  cycles. 

Computer  programs  using  Hardin's  soil  model  employ  an  iterative 
scheme  to  represent  the  non-linear  soil  behavior.  During  each  iter¬ 
ation,  a  stiffness  matrix  is  formed  and  nodal  displacements  are  com¬ 
puted.  The  nodal  displacements  are  used  to  calculate  element  stresses 
and  strains.  At  the  end  of  each  iteration  the  shear  modulus  of  the  soil 
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Is  modified  Co  provide  an  Improved  agreement  between  element  stresses 
and  strains.  Between  three  and  five  iterations  are  usually  sufficient 
to  provide  an  accurate  analysis  using  Hardin's  model.  The  stiffness  of 
the  soil  decreases  with  successive  iterations. 

Other  relations  established  included  a  non-dimensional  relationship 
for  Che  unloading  shear  modulus,  G^,  and  the  unloading  shear  strain,  y  . 
For  a  silty  sand  and  a  silty  clay  soil  the  permanent  set  strain,  S,  was 
related  to  the  number  of  load  applications  N,  by 

S  -  1  +  log^N  (4) 

where  y/y  was  less  than  0.3. 
r 

Y.  T.  Chou  and  R.  H.  Labbetter  (1973)  reported  on  the  behavior  of 
flexible  airport  pavements  underloads.  He  reported  about  theory  and 
experiment.  One  purpose  of  their  study  was  to  reduce  and  analyze  the 
deflection,  stress  and  strain  data  obtained  from  a  series  of  multiwheel 
aircraft  flexible  pavement  tests.  The  second  purpose  was  to  establish 
relations  between  load  and  pavement  response. 

It  was  reported  that,  since  pavement  materials  possess  non-linear 
material  characteristics,  the  use  of  non-linear  finite  element  solutions 
seemed  to  be  an  obvious  solution  scheme.  However,  the  development  of 
the  constitutive  relations  at  Chat  time  could  not  adequately  express  the 
complex  stress  changes  in  a  pavement  system.  Despite  the  nonlinearity 
of  pavement  materials,  the  relations  between  stress  and  displacements 
were  linear  in  some  aspects  and  the  principal  of  superposition  was 
reasonably  valid. 

Brown  (1974)  did  some  research  into  repeated  loading  testing  of  a 
granular  material  in  a  triaxial  chamber.  He  established  that,  for  a 
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dynamic  deviator  stress,  q,  equal  to  twice  the  mean  value,  p,  both  the 

4 

resilient  strains  and  permanent  strains  changed  very  little  after  10 
cycles.  The  permanent  strain,  e^,  was  linearly  related  to  the  applied 
stress  by  ■  0.01  q/a^,  where  is  confining  pressure.  The  resilient 

strain,  ,  was  related  to  the  applied  stress  by 


(p+.25c3J 


microstrain 


McLean  and  Brooker  (1975)  did  a  study  to  estimate  the  permanent 
deformations  in  asphaltic  concrete  due  to  repeated  traffic  loadings.  In 
their  article,  they  discussed  a  methodology  to  estimate  permanent 
deformations  (rutting)  in  pavements  from  repeated  load  triaxial  com¬ 
pression  tests  and  creep  tests.  Analytical  calculations  were  made  to 
predict  rut  development  using  the  elastic  layer  program  CHEV5L.  A 
viscoelastic  analysis  of  creep  prediction  was  done  using  the  VESYSII 
program. 

In  the  laboratory  program,  repeated  triaxial  load  tests  were  made 
on  asphaltic  bound  materials  over  a  range  of  temperatures.  The  tech¬ 
nique  for  estimating  rut  depth  involved  fitting  a  curve  through  the 
results  of  the  repeated  load  tests.  A  third  order  polynomial  was 
evolved  relating  the  plastic  strains  (corrected  for  volume  change),  e^, 
to  the  number  of  stress  applications,  N,  by 

Ep  -  CQ  +  C1log(N)  +  C2(logN)2  +  C3vlogN)3.  (6) 

C^,  C2  and  C3  were  experimental  constants  with  mean  values  of  0.85, 
0.013,  and  -0.14  respectively.  C_  was  the  logarithm  of  the  permanent 
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strains  after  the  initial  load  application,  c^.  was  influenced  by 
temperature,  stress,  and  strain  levels.  The  following  relationship 
was  developed  to  show  the  influence: 


e 


P 

1 


K<VV  • 


(7) 


where  ad  is  the  difference  in  axial  and  radial  stress  in  triaxial 
compression,  is  elastic  strain,  and  K  and  n  were  experimental  con¬ 
stants  dependent  of  temperature.  Axial  creep  compliance  was  established 
through  tension  tests  at  various  temperatures.  The  creep  compliance 
ifi£(t)  as  a  function  on  temperature,  t,  was  defined  as: 


♦e<t) 


eft) 

a 


(8) 


where  e  (t)  is  the  axial  strain,  and  o,  is  defined  above, 
a  a 

The  analytical  procedure  used  an  elastic  analysis  of  a  two  layer 
system,  asphalt  over  subgrade.  The  wheel  load  was  applied  over  a 
circular  area.  The  strains  in  the  pavement  system  were  related  to  the 
triaxial  tests  by  assuming  a,j  “  0X  ~  where  is  the  vertical  stress 
and  is  the  horizontal  stress.  The  permanent  deflections,  6^,  were 
computed  by  dividing  the  pavement  into  a  number  of  layers  of  thickness 
AZ^.  The  permanent  strains  in  each  layer  were  computed  using  the 
methodology  described  above.  The  sum  of  the  permanent  deflections, 
x  AZ^  of  each  layer  were  summed  to  yield  the  total  surface 
deflections,  or 

No.  of  layers 
6P  -  Z  (eP'AZ.) 

i-1  1  1 


(9) 
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Poulsen  (1978)  reported  on  the  laboratory  testing  of  cohesive 
subgrades.  Potential  modes  of  distress  for  subgrates  were  considered  to 
be  directly  or  indirectly  related  to  (1)  elastic  or  resilient  loadings 
caused  by  traffic  loadings,  (2)  permanent  or  plastic  deformations,  and 
(3)  plastic  deformations  caused  by  climatic  or  environmental  changes,  as 
changes  in  soil  suction,  freeze-thaw,  or  differential  settlement. 

Poulsen  directed  his  research  to  determining  the  relation  between 
plastic  deformations  and  traffic  loadings. 

His  procedure  was  to  take  undisturbed  samples  from  6  countries 
including  England,  West  Germany  and  the  United  States.  The  samples  were 
subjected  to  triaxial  repeated  load  tests.  A  constant  confining  pres¬ 
sure  and  a  vertical  static  preload  were  applied  to  all  samples.  The 
minor  principal  stress  ranged  from  5  to  10  kPa  (0. 7-1.4  psi) .  The  major 
principal  stress  was  approximately  twice  the  minor  principal  stress.  A 
vertical  dynamic  stress,  was  superimposed  on  the  static  stress 

state. 


To  establish  a  true  picture  of  the  process  of  permanent  strain 


accumulations  for  a  given  a  large  number  of  cyclic  load  tests  were 

done.  A  basic  simple  expression 


e 

P 


(10) 


was  presented  to  represent  the  test  results,  e. ,  e.,  and  e,  were 
material  constants.  N  was  the  number  of  load  applications  in  a  single 
stage  test.  8  was  the  ratio  °{jyn/a(jyT1  f  called  the  degree  of  failure, 
where  ^  was  the  dynamic  deviator  stress  at  triaxial  failure  for  N  - 

100,000  cycles.  The  equation  was  considered  valid  for  N  greater  than 
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100  cycles  and  8  less  chan  one.  The  proposed  equation  had  no  apparent 
correlation  with  general  soil  parameters. 

Majidzadeh  et  al.  (1978)  reported  on  the  rutting  evaluation  of 
subgrade  soils  in  Ohio.  The  object  of  the  report  was  to  study  rutting 
experimentally  by  using  laboratory  prepared  samples  subjected  to 
uniaxial  dynamic  testing.  The  following  rutting  model  was  appropriate 
to  the  five  types  of  silty  soils  used  in  the  study: 

-e  /N  -a/  (11) 


where  is  the  permanent  strain,  N  is  the  number  of  stress  appli¬ 
cations,  m  is  the  slope  of  the  straight  line  relation  between  log(tp/N) 
and  logN,  A  is  at  N  equal  to  1.  The  parameter  m  varies  between  0.85 
and  .90. 


The  rutting  parameter  A  was  established  as  a  function  of  E*,  the 
dynamic  modulus,  by 


A  *  R(E*)-C  exp  {aapl/oult> .  (12) 

R  and  C  are  experimental  constants.  o  ^  is  the  applied  stress  and 
is  the  ultimate  stress  from  an  unconfined  strength  test. 

Frazier,  Parker,  et  al.  (1979)  reported  a  structural  design  pro¬ 
cedure  for  rigid  airport  pavements.  They  selected  as  a  response  model 
for  the  pavement  system  the  computer  code  BISTRO.  The  code  models  the 
pavement  as  layered  elastic  materials.  The  design  criteria  was  based  on 
the  tensile  stress  of  concrete  as  computed  by  the  program  and  the 
strength  of  concrete  as  measured  by  flexural  beam  tests.  The  limiting 
stress  criteria  was  related  to  the  number  of  coverages  through  fatigue 


testing  and  converting  the  traffic,  load  of  multiple  types  of  aircraft  to 
an  equivalent  loading  of  a  single  aircraft. 

They  also  reported  in  the  Appendix  the  engineering  behavior  of 
pavement  materials.  They  reported  on  many  different  materials  in 
pavements  and  their  general  engineering  behavior  due  to  repeated 
loadings.  A  summary  is  given  below. 

The  effects  of  repeated  loadings  for  concrete  were  generally 
handled  indirectly  by  fatigue  relationships.  However,  the  magnitude  of 
stress  sustained  before  cracking  was  a  function  of  the  number  of  repe¬ 
titions,  and  the  magnitude  of  this  stress  decreased  as  the  number  of 
repetitions  increased. 

Bituminous  materials  were  reported  to  be  viscous  and  temperature 
dependent.  Their  behavior  was  affected  to  a  greater  extent  by  the  rate 
of  loading  than  by  chemical  composition.  Complete  characterization 
requires  different  rates  of  loading  so  runways,  taxiways  and  aprons 
could  be  adequately  characterized.  From  a  practical  point  of  view,  the 
rate  of  loading  and  temperature  variations  are  limited.  It  is  possible 
to  include  the  effects  of  the  loading  rate  and  temperature  variations 
during  repeated  load  tests. 

Granular  bases  were  reported  to  be  extremely  difficult  to  char¬ 
acterize.  The  state  of  stress,  particularly  the  confining  pressure,  was 
the  dominating  factor  in  determining  the  load  deformation  properties. 
Repeated  load  properties  were  different  from  static  properties.  Re¬ 
peated  loads  tend  to  increase  the  stiffness,  provided  a  progressive 
shear  failure  does  not  occur.  Poisson's  ratio,  v,  was  reported  to  vary 
during  shear  failure  and  the  direction  of  change  was  dependent  on  the 
initial  density.  For  low  density  soils,  densif ication  occurs  during 
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shear,  and  v  decreased.  Poisson's  ratio  would  increase  if  the  material 
increased  in  volume  during  failure.  In  general,  v  reaches  a  relative 
constant  value  during  repeated  load  tests  and  will  remain  constant 
unless  shear  failure  occurs. 

Subgrades  were  reported  to  be  the  components  of  a  rigid  pavement 
system  most  affected  by  repeated  load  applications.  Host  natural  soils, 
when  in  the  saturated  condition,  behave  primarily  like  cohesive  soils. 

Snaith  et  al.  (1980)  reported  on  flexural  pavement  analysis.  Their 
research  was  applied  to  vehicular  pavements.  They  used  the  finite 
element  program  DEFPAV  as  a  pavement  model  and  were  able  to  simulate 
creep  behavior,  crack  propagation  through  the  bitumunious  layer,  and 
compute  the  transient  deflections  with  different  cracking  ratios.  They 
also  calculated  the  growth  of  rut  depths  and  compared  their  calculations 
to  a  test  section. 

The  creep  characterization  of  the  pavement  materials  was  derived 
from  repeated  load  triaxial  tests.  The  creep  equations  for  the  subgrade 
soils  used  a  non-iterative  procedure  to  relate  the  induced  permanent 
strain  to  the  imposed  stress.  The  surface  deflections  were  computed  by 
summing  the  displacements  in  the  elements  by 
No.  of  elements 

6  -  Z  eh,  (13) 

i-1  1  1 

where  was  the  element  permanent  strain  and  h^  was  the  height  of  the 
element.  Their  work  was  extended  to  overlay  pavement  design. 

Selig,  et  al.  (1981)  developed  a  methodology  for  computing  the 
cumulative  permanent  deformations  of  a  railroad  track  structure  due  to 
traffic  loading.  They  wanted  to  use  this  procedure  as  a  basis  for 
determining  the  maintenance  life.  They  used  a  computer  model  called 
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GEOTRACK.  Available  Information  on  the  elastic  and  inelastic  properties 
of  ballast  and  subgrade  materials  was  reviewed.  Static  and  cyclic 
triaxlal  tests  were  carried  out  on  granite  ballast.  Their  methodology 
of  computing  permanent  strain  was  compared  to  permanent  strains  on  the 
FAST  track. 

They  concluded  that  the  vertical  permanent  strains  after  the  first 
loading  cycles  could  be  predicted  by  two  methods.  One  was  to  use  the 
results  of  cyclic  test  results  which  depict  the  permanent  strains  as  a 
function  of  the  ratio  Aq/q^  where 


V°3 


V°3 


qf  (a^-a^Jf  2c  cos4>  +  2a2siI*<|>  v  ' 

1  -  sini)> 

is  the  triaxlal  normal  pressure,  is  the  confining  pressure,  c  is 
the  cohesion,  $  is  the  angle  of  internal  friction,  and  the  subscript  f 
refers  to  failure  conditions.  Typical  lab  curves  are  shown  in  Figure  6. 
The  induced  strain  is  a  function  of  the  soil  properties  c  and  <(>,  and  the 
stress  level. 

The  second  procedure  for  computing  the  permanent  strain  after  one 
cycle,  e^,  was  to  use  hyperbolic  parameters  from  static  triaxlal  tests. 
The  axial  strain  predicted  by  hyperbolic  modeling,  ev>  can  be  computed 


°1  °3/El 

R  (l-sin$) 

2c  cos$  +  2a.jSin$ 


where  R^  was  defined  as  the  failure  ratio  and  was  the  initial  tangent 

N 

modulus  equal  to  Ka^,  where  K  and  N  were  experimentally  determined  con¬ 
stants.  The  value  of  can  be  computed  by  Equation  16  from  if  the 
resilient  strain  after  the  first  cycle,  e^,  is  known. 
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(16) 


e  has  to  be  determined  from  testing  and  has  shown  to  be  related  to 
Aq/q^  as  shown  in  Figure  7. 

The  value  of  the  permanent  strain  after  N  cycles,  e^,  was  related 
to  the  first  cycle  permanent  strain  by  a  laboratory  analysis  of  cyclic 
triaxial  testing.  Both  a  linear  and  non-linear  equation  were  presented 
for  ballast. 


(1+  19  log1Q(Nl)  linear  equation  (17a) 

ew  - 

(.85+. 38  logN+e^( .05)-. 08  logN)  non-linear  equation  (17b) 

A  comparison  of  the  computed  strain  values  and  actual  strains  in  FAST 
track  compared  favorably. 


Conclusion 

A  brief  description  of  the  methods  of  computing  residual  displace¬ 
ments  and  the  knowledge  of  the  behavior  of  pavement  materials  as  deve¬ 
loped  over  the  past  twenty  years  has  been  presented.  The  current  method 
of  computing  residual  displacements  involves  the  use  of  finite  element 
or  elastic  layer  theory  computer  programs  to  compute  the  initial  set 
strains  due  to  the  first  loading.  Final  set  strains  are  related  to  the 
initial  set  strains  via  actual  testing  of  the  material  under  triaxial 
conditions.  The  final  rut  depths  are  computed  using  Equation  2.  These 
types  of  procedures  are  valid  as  long  as  the  pavement  materials  do  not 
undergo  shear  failure.  Monograph  solutions  are  also  available  but  were 


not  mentioned  in  the  Literature  Review. 
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CHAPTER  III 

Methodology  for  Computing  Residual  Displacements 

Introduction 

In  this  chapter  two  methodologies  for  computing  residual  displace¬ 
ments  in  surfaced  and  unsurfaced  pavement  systems  are  presented.  All 
these  methods  involve  the  use  of  finite  element  computer  codes.  Trial 
solutions  for  two  of  the  methods  were  attempted  with  varying  degrees  of 
success.  To  develop  each  procedure  to  a  state  where  it  could  be  used  as 
an  everyday  tool  was  beyond  the  scope  of  this  project.  Sample  solutions 
were  attempted  to  determine  the  feasibility  of  each  procedure  rather 
than  to  fully  develop  an  excellent  software  package  for  each  method. 

Procedure  No.  1 

The  procedure  discussed  in  this  section  closely  approximates  the 
state  of  art  of  residual  displacement  computation  as  revealed  by  the 
literature  review.  It  has  three  basic  steps.  First,  finite  element 
computer  codes  were  used  to  predict  the  permanent  strains  induced  by  one 
pass  of  an  aircraft  tire  over  a  pavement  section.  Secondly,  the  verti¬ 
cal  strains  induced  by  repeat  load  applications  was  estimated  by  using 
laboratory  triaxial  repeat  load  data.  Thirdly,  the  surface  deflections 
were  computed  by  summing  the  change  of  height  of  each  element.  Mathe¬ 
matically  this  was  expressed  by 

Num 

S  -  Z  z  *h  (18) 

i-1 

where  S  is  the  surface  deflection,  e  is  the  element  strain  induced  by 
repeated  loads,  h^  is  the  height  of  each  element,  and  Num  is  the  number 
of  elements  beneath  the  surface  where  the  rut  depth  is  being  calculated. 


Two  different  computer  programs  were  used  to  compute  the  initial 
residual  strains  after  one  wheel  load.  These  were  the  Bomb  Damage 
Repair  (BDR)  code  (Rude,  1980)  and  the  Culvert  Analysis  and  Design 
(CANDE)  code  (Katona,  1976).  These  codes  were  selected  because  they 
contained  the  non-linear  material  law  developed  by  Hardin  (1971). 

The  BDR  code  can  represent  pavement  systems  either  with  an  axisym- 
metric  mesh  or  with  a  mesh  composed  of  prismatic  solid  elements.  Only 
the  axisymmetrical  mesh  was  used.  The  CANDE  program  was  developed  to 
design  and  analyse  buried  culverts  using  two  dimensional  plane  strain 
elements  for  soil  and  beam  column  elements  for  the  culvert.  Pavement 
systems  were  represented  using  the  quadlateral  elements  and  omitting  the 
beam  column  elements. 

The  calculations  using  the  BDR  code  were  done  at  Tyndall  Air  Force 
Base  during  the  summer  of  1980.  The  program  was  written  for  CDC 
machines.  A  great  deal  of  time  was  used  to  attempt  to  adopt  the  BDR 
code  to  the  host  IBM  computer  at  Virginia  Tech  (VPI&SU) .  The  program 
was  not  successfully  running  towards  the  close  of  the  project,  so  much 
of  the  sample  problems  were  done  using  CANDE.  CANDE  was  the  only  other 
program  known  to  contain  Hardin's  soil  model. 

According  to  the  literature  review,  the  greatest  permanent  strains 
occur  during  the  initial  application  of  the  wheel  load.  Most  pavement 
materials  exhibit  a  non-linear  relation  between  stress  and  strain.  A 
non-linear  stress-strain  law  seemed  the  most  appropriate  type  of  law  to 
use  to  represent  the  initial  load.  Hardin's  law  was  chosen  to  represent 
the  pavement  materials  because  of  its  non-linear  nature  and  correlation 
to  traditional  material  characterization  parameters  as  void  ratio 
density  and  plasticity  index. 
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A  typical  stress  strain  curve  for  a  non-linear  material  undergoing 

cyclic  loading  is  shown  in  Figure  2.  Shown  in  the  figure  are  the 

initial  tangent  modulus  E^,  the  unloading  modulus,  E^,  the  permanent 

strain  after  one  cycle,  e, ,  and  the  permanent  strain  after  N  cycles,  S„. 

i  N 

The  literature  review  has  shown  that  S„  can  be  defined  in  terms  of  e,  by 

N  1 

a  logrithmic  expression  or  a  hyperbolic  power  expression. 

For  the  equation  to  be  used  to  compute  the  permanent  set  strain, 
the  initial  set  strain  e^  must  be  determined.  Using  a  finite  element 
computer,  this  would  involve  loading  the  mesh  using  a  variable  modulus 
and  then  removing  the  effects  of  the  load  using  the  unload  modulus. 

None  of  the  available  computer  codes  contained  the  programming  to  apply 
and  successfully  remove  a  load. 

The  following  method  was  adopted  to  compute  e^  instead.  The 
aircraft  wheel  load  was  applied  to  the  pavement  using  non-linear  consti¬ 
tutive  material  models.  This  analysis  corresponds  to  the  loading  phase. 
To  approximate  the  unloading  phase,  the  same  wheel  loading  was  applied 
to  the  same  grid,  but  all  materials  were  represented  by  the  initial 
tangent  modulus.  The  initial  tangent  modulus  was  assumed  to  be  equal  to 
the  unload  modulus.  In  reality  this  is  only  approximately  correct.  The 
element  strains  in  the  unloading  phase  were  subtracted  from  the  element 
strains  obtained  during  the  loading  phase  to  yield  the  permanent  set 
strains  after  one  cycle. 

Several  attempts  were  made  to  simulate  unloading  by  using  two 
loading  increments  with  the  CANDE  program  with  little  success.  The  load 
was  applied  during  the  first  loading  increment  and  a  load  of  the  same 
magnitude,  but  opposite  sign,  was  applied  during  the  second  loading 
increment  to  zero  out  the  first  load.  Unrealistic  strains  and 
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displacements  were  obtained.  Hardin's  lawwas  used  for  these  calcu¬ 
lations.  The  CANDE  program  reduced  the  stiffness  of  any  material  when 
tensile  stresses  were  applied.  This  stiffness  reduction  produced  the 
unrealistic  results.  This  is  why  only  compressive  forces  w«.re  applied 
when  using  CANDE. 

Other  limitations  were  also  discovered  using  the  CANDE  program  to 
compute  the  first  cycle  residual  strain  in  clay  soils.  Hardin's  law  is 
programmed  to  perform  one  to  five  Iterations  to  find  the  appropriate 
secant  shear  modulus  for  a  particular  state  of  stress.  The  first  cycle 
corresponds  to  using  the  maximum  shear  modulus.  An  analysis  of  "load" 
and  "unload"  calculations  showed  that  the  strains  for  the  "unload"  case 
were  greater  than  the  "load"  case.  This  seemed  theoretically  impos¬ 
sible,  so  fine  grain  soils  were  modeled  as  linear  materials.  For  the 
unload  case,  an  initial  tangent  Young's  modulus  was  used,  and  for  the 
load  case  an  appropriate  secant  Young's  modulus  was  used. 

Comparison  of  Procedure  with  "Buckshot  Clays"  Tests 

G.  W.  Tumage,  et  al.  (1974)  perform  large  scale  laboratory  tests 
to  correlate  the  ruts  produced  from  two  different  types  of  aircraft 
tires  and  three  standard  military  trucks.  These  tests  were  done  on  the 
"Buckshot  clay"  at  various  compacted  densities  that  were  placed  in 
indoor,  concrete-lined  pits  6  feet  deep  and  11.7  feet  wide  by  13  feet 
long.  A  high  strength  clay  was  placed  within  24  inches  of  the  top  of 
the  pit  and  enough  soil  was  compacted  at  specified  densities  and 
moisture  contents  using  a  40,000  lb  self-propelled,  multiwheel  roller  to 
fill  the  rest  of  the  tank. 


mmm wm*  ■*»-*ti*m**v:  ^*^*w****** 
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The  buckshot  clay  is  a  Mississippi  River  alluvium.  It  is  a  highly 
plastic,  almost  purely  cohesive  soil.  It  has  a  CH  Unified  Soil  Classi¬ 
fication.  For  the  comparison  study,  test  data  with  the  soil  compacted 
at  a  dry  density  of  104.2  pcf,  at  22.5  percent  moisture  content,  and  a 
degree  of  saturation  of  98.6  percent  was  used.  The  soil  was  trafficked 
with  a  C130  aircraft  tire  (20-20,  22  pr)  with  a  total  load  of  approxi¬ 
mately  35384  lbs  and  a  tire  pressure  of  100  psi.  The  aircraft  tire  was 
mounted  on  a  carriage  that  spanned  the  tank.  The  carriage  restricted 
the  tire  loading  to  a  single  path  through  the  soil.  An  average  force  of 
1549  lbs  was  used  to  tow  the  tire  through  the  soil.  The  results  of  the 
rut  depth  versus  number  of  repetitions  is  shown  in  Figure  8.  The  slope 
of  a  curve  fitted  through  Turnage's  data  is  0.546. 

For  the  computer  modeling,  the  initial  loading  and  unloading  of  the 
tire  were  simulated  by  applying  a  100  psi  pressure  to  the  surface  of  the 
finite  element  grid  shown  in  Figure  9.  The  CANDE  program  was  used  and 
the  soil  was  modeled  as  a  linear  elastic  material  with  a  Young's  modulus 
of  “  5700  psi  for  loading  and  E^  ■  11400  psi  for  unloading. 

Poisson's  ratio  of  the  soil  was  0.43.  The  computed  permanent  strains 
after  one  cycle  were  obtained  by  subtracting  the  respective  element 
strains  computed  using  E  from  the  element  strains  computed  using  E  . 

U  la 

The  total  rut  Jepth  was  computed  using  Equation  18. 

Turnage's  report  did  not  contain  repeated  load  test  data  on  the 
''Buckshot  Clay".  They  were  obtained  by  assuming  an  empirical  relation 
between  accumulated  strains  and  the  initial  permanent  strains.  The 
assumption  was  made  that  the  permanent  strains  would  accumulate  as  the 
same  rate  as  the  total  rut  depth.  Using  the  data  provided  in  Figure  8 
and  the  mathematical  formulation  suggested  by  Poulsen,  the  following 
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expression  was  obtained:  -  e^N  ,  where  B  -  0.546,  N  ■  number  of  load 

applications,  -  permanent  strains  after  one  cycle  and  e^  -  accumu¬ 
lated  set  strain. 

Table  1  summarizes  the  calculations  for  the  finite  elements  beneath 

the  center  of  the  wheel  pressure.  The  initial  rut  depth  of  0.14  in 

agree  with  the  extrapolated  test  value  of  0.155  in.  The  final  rut  depth 

of  1.73  in  agree  well  with  the  actual  value  of  1.77  in  after  100  passes. 

Additional  data  from  Tumage's  tests  were  plotted  in  a  graph  of  rut 

depth  versus  number  of  applications.  The  test  results  corresponded  to 

the  application  of  a  C-130  tire  with  the  "Buckshot  Clay"  compared  at 

different  densities  at  different  moisture  contents.  See  Figure  10.  The 

slopes  of  each  curve  are  approximately  equal. 

Tumage  also  reported  the  results  of  a  cone  penetrometer  test  for 

each  filling  of  the  testing  tank.  Figure  10  also  shows  the  results  of 

plotting  the  rut  depth  after  the  first  pass  versus  the  cone  index.  The 

cone  index.  Cl,  was  the  force  (in  lbs)  per  unit  area  of  the  cone  base 

required  to  penetrate  the  soil  in  a  vertical  direction  at  a  rate  of  72 

in/min.  The  cone  was  a  right  angle  cone  with  a  30°  apex  angle  and  a  0.5 
2 

in  base  area.  These  results  would  tend  to  indicate  that  a  cone  pene¬ 
trometer  test  could  be  used  to  estimate  the  first  rut  depth,  h^,  and  the 

g 

expression  h  -  h^N  could  be  used  to  estimate  the  rut  depth,  h,  produced 
after  N  loadings.  If  the  hypothesis  is  carried  further,  Che  same 
relations  could  hold  for  estimating  rut  predictions  via  computer 
analysis. 


ilTWiifiw  »i  i  •**••*-**■*“’* 
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Evaluation  of  Crushed  Limestone  Crater  Repair 

The  current  method  for  an  expedient  repair  of  a  runway  crater 
involves  four  steps.  These  are:  (1)  pushing  debris  into  the  crater  up 
to  24  inches  below  the  pavement  surface,  (2)  placing  select  fill 
(crushed  limestone)  into  the  crater  and  overfilling  by  approximately  6 
inches,  '3)  compaction  using  a  self  propelled  vibratory  roller  and  (4) 
placement  of  a  Foreign  Object  Damage  (FOD)  cover.  The  Engineering  and 
Services  Laboratory  at  Tyndall  Air  Force  Base,  Florida  have  conducted 
several  simulated  repairs  and  have  trafficked  them  with  an  F-4  load 
cart.  One  repair  was  briefly  reported  by  Rude  (1980)  and  three  were 
reported  in  detail  by  Knox  (1980). 

The  repair  reported  by  Rude  (i980)  was  made  by  placing  24  in  of 
crushed  limestone  into  a  simulated  crater.  At  the  bottom  of  the  crater 
was  a  clay  soil.  The  repair  was  trafficked  by  1440  coverage  of  an  F-4 
load  cart  distributed  among  twelve  paths.  The  maximum  number  of  passes 
on  any  strip  was  150  passes  of  the  load  cart.  Figure  11  shows  a  simu¬ 
lated  picture  of  the  crater  profile  after  1440  coverages.  The  profile 
was  obtained  by  subtracting  the  elevation  at  selected  points  in  the 
repair  after  trafficking  from  the  elevation  before  trafficking.  The 
repair  was  assumed  to  be  level  after  the  compaction  process  was  com¬ 
pleted.  The  average  rut  depth  for  the  area  that  received  150  passes  was 
1.87  in. 

Rude  (1980)  used  the  Bomb  Damage  Repair  (BDR)  Computer  Program  to 
estimate  the  ruts  produced  by  an  F-4  wheel  loading.  THe  BDR  program 
used  Hardin's  law  to  represent  the  stress-strain  behavior  of  both 
materials.  Using  the  procedure  discussed  in  the  theory  section,  the 
permanent  strains  after  one  pass  of  an  F-4  load  cart  as  calculated  from 
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Che  BDR  calculations  was  1.12  In.  The  relationship  that  the  accumulated 
set  strains  ■  e^(l+.19  logN),  where  -  first  cycle  set  strains  and 
N  -  number  of  load  applications,  was  used  to  predict  the  rut  depth  after 
150  passes.  The  calculated  rut  depth  was  1.61  in.  This  prediction  was 
13.5%  of  the  average  value.  Only  the  strains  in  the  finite  elements 
within  the  top  12  in.  of  the  select  fill  and  the  top  10  in.  of  the  clay 
had  any  significant  contribution  to  the  rut  development.  Significant 
element  strains  produced  a  permanent  settlement  of  0.001  in.  or  greater. 

The  crater  repairs  reported  by  Knox  (1980)  were  conducted  for  the 
purpose  of  field  testing  the  rapid  repair  technique  and  to  evaluate  the 
repairs  with  a  F-4  load  cart.  Three  crater  repairs  were  properly  made 
using  crushed  limestone  placed  on  debris  backfill.  Crater  repair  number 
three  had  an  average  of  2  1/2  ft  of  crushed  limestone,  crater  repair 
number  five  had  approximately  36  in.  of  select  fill,  and  crater  number 
six  had  between  18  in.  to  30  in.  of  select  fill.  Each  repair  received 
1440  coverages  of  an  F-4  load  cart.  Crater  3  was  repaired  after  100 
coverages  by  additional  amounts  of  limestone  due  to  consolidation  of  the 
repair.  Repairs  number  five  and  six  received  96  coverages  of  the  F-4 
load  cart  immediately  after  being  made.  An  additional  1344  coverages 
were  made  during  the  next  16  days.  The  delay  was  due  to  repairs  re¬ 
quired  to  the  F-4  load  cart.  The  maximum  rut  depth  due  to  repairs 
three,  five  and  six  were  1.7  in.,  0.8  in.,  and  1.4  in.,  respectively. 
Knox  reported  that  the  time  delay  in  trafficking  the  repairs  probably 
resulted  in  the  abnormally  small  surface  deflections  for  repairs  five 
and  six. 

The  CANOE  computer  program  was  used  to  simulate  crater  number  six. 
The  depth  of  the  select  fill  was  assumed  to  be  24  inches.  The  computer 
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program  was  used  Co  calcualte  the  element  strains  due  to  applying  the 
F-4  wheel  load  and  removing  the  load.  The  difference  In  these  strains 
was  the  permanent  strains  induced  after  the  first  cycle.  Hardin's  soil 
model  was  used  to  represent  the  gravel.  Five  iterations  were  used  for 
the  loading  and  one  iternation  was  used  for  the  unloading.  This  push- 
back  was  modeled  as  a  linear  material. 

Knox  did  not  report  field  or  laboratory  information  with  regard  to 
the  "compacted"  pushback  material.  Green  (1978)  summarizes  elastic 
constants  for  airport  material.  Green  reports  on  work  by  Foster  and 
Heukelom  that  says  the  subgrade  modulus  for  soil  can  vary  from  1130  psi 
for  peats  to  127,000  psi  for  moraine  deposits.  No  information  was 
listed  regarding  the  physical  properties  of  pushback  compacted  according 
to  the  BDR  repair  procedure.  Without  actual  values  for  the  pushback 
material  only  ball-park  estimates  of  rutting  behavior  were  determined. 

A  wide  range  of  values  for  the  elastic  modulus  was  chosen  to 
represent  the  behavior  of  the  pushback  material.  The  results  of  the 
series  of  calculations  are  shown  in  Figure  12.  Plotted  on  the  vertical 
axis  are  the  maximum  surface  vertical  deflections  (in  inches)  and 
plotted  on  the  horizontal  axis  is  Young's  modulus  of  the  pushback  (psi) 
on  a  log  scale.  The  upper  curve  represents  the  surface  deflections  for 
the  loading  case  and  the  lower  curve  represents  the  unloading  case.  To 
compute  the  initial  residual  displacement  after  one  cycle,  the  surface 
displacements  according  to  curve  B  were  subtracted  from  the  surface 
displacement  of  curve  A.  These  results  agreed  well  when  equation  18  was 
used  to  compute  the  first  cycle  permanent  deflections. 

Data  analysis  indicates  that  in  the  limiting  case  where  the  modulus 
of  the  pushback  is  exceptionally  high,  the  select  fill  deflects  1.30  in. 
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after  one  loading,  and  has  a  permanent  rut  deflection  of  1.85  inches 

after  150  passes.  If  the  pushback  was  highly  compacted  with  an  unload 

modulus  of  60,000  and  a  load  modulus  of  6,000  psi,  the  first  cycle 

permanent  deflection  of  the  repair  would  be  2.37  in.  and  the  maximum  rut 

depth  after  150  cycles  would  be  3.72  in.  In  this  estimate,  the  select 

fill  was  assumed  to  behave  according  to  the  relationship  -  Cp^(l-t-.19 

logN)  and  the  pushback  according  to  the  relationship  ■  e^d+logN) 

where  e  ,  is  the  residual  element  strain  after  one  cycle  and  e  is  the 
pi  P 

residual  strain  after  N  cycles.  The  permanent  surface  delections  were 
computed  accroding  to  equation  18. 

The  computed  results  were  much  larger  than  reported  by  Knox.  The 
average  depth  of  the  select  fill  could  have  been  less  than  the  24  in. 
assumed.  The  plane  strain  elements  in  the  CANDE  program  may  not  be 
applicable  for  this  crater  repair.  Improved  calculations  might  result 
if  a  finer  finite  element  mesh  has  been  used.  In  general,  however,  an 
estimated  range  of  deflections  from  1.3  in.  to  3.72  in.  does  not  seem  to 
be  unreasonable. 

A  possible  field  check  to  determine  proper  compaction  of  pushback 
material  would  be  to  employ  a  cone  penetrometer.  Sanglerat  (1972)  has 
reported  a  number  of  research  projects  where  correlations  have  been 
established  and  confirmed  by  other  researchers  for  compaction  control 
using  the  one  penetrometer.  Computer  analysis  of  a  particular  repair 
plan  could  determine  limiting  values  for  subgrade  modulus  and  densities 
for  the  pushback.  Correlations  could  be  established  between  different 
densities  of  pushback  and  cone  penetration  values.  The  cone  penetro¬ 
meter  could  then  be  used  as  a  field  check  to  determine  if  sufficient 
compaction  had  been  achieved. 
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Applications  to  Airport  Pavement 

A  sample  problem  was  generated  to  apply  the  method  for  computing 
residual  displacements  to  a  concrete  runway.  The  pavement  section 
consisted  of  a  concrete  over  a  granular  base  with  a  stiff  clay  as  a 
subbase.  The  properties  and  the  height  of  each  layer  as  idealized  in 
the  finite  element  grid  are  given  in  Table  2. 

A  series  of  calculations  were  made  to  determine  the  maximum  rut 
depth  under  an  F-4  wheel  loading.  The  process  previously  described  was 
used.  The  results  of  these  calculations  are  shown  in  Figure  13.  The 
horizontal  axis  is  the  number  of  repetitions,  and  the  vertical  axis  is 
the  maximum  rut  depth.  The  concrete  was  c’.sumed  to  settle  along  with 
the  pavement  base  and  subgrade.  The  rate  of  induced  strain  due  to  load 
repetitions  for  each  material  is  shown  in  Table  3. 

In  an  actual  runway,  aircraft  do  not  land  or  take  off  along  a 
definite  path.  According  to  HoSang  (1978),  the  aircraft  traffic  distri¬ 
bution  can  be  represented  by  a  theoretial  normal  distribution.  The  mean 
value  is  slightly  offset  from  the  centerline.  A  typical  histogram  for 
the  lateral  distributions  of  aircraft  traffic  during  takeoffs  is  shown 
in  Figure  14. 

To  illustrate  the  process  of  computing  the  residual  displacements, 
a  series  of  F-4  wheel  loadings  for  a  given  number  of  total  applications 
were  distributed  across  the  pavement  section  according  to  the  distri¬ 
bution  shown  in  Figure  14.  The  accumulated  deflections  were  then 
calculated  for  the  prescribed  number  of  load  applications.  The  total 
residual  displacements  for  a  total  of  200,  20,000,  and  200,000  load 
applications  are  presented  in  Table  4  and  illustrated  in  Figure  15. 
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Table  3 

Rate  of  Induced  Strain  Due  to  Load  Repetitions 


Layer 


Composition 


Strain  Accumulation  Equation 


1 

Concrete 

None 

2 

Gravel 

e  -  e,  (1  +  .19  logN) 
n  l 

3 

Clay 

„0.55 
e  -  e.N 
n  1 

I  ■  number  of 

load  applications 

Percent  Distribution 


Distance  (ft)  from  Runway  Centerline 


Figure  14.  Typical  Histogram  for  Lateral  Distribution  of 

Aircraft  Traffic  During  Takeoffs,  after  HoSang  (1978) 


Permanete  Deflections  Across  Concrete  Runway 
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The  illustrative  example  is  very  simplistic  in  that  it  represents 
that  of  only  the  application  of  one  wheel  loading  and  not  the  total 
wheel  loading  of  an  F-4  aircraft  with  its  own  landing  gear  configur¬ 
ation.  An  algorithm  could  easily  be  written  to  distribute  the  actual 
number  of  wheel  loadings  to  a  particular  location  on  the  pavement.  If  a 
particular  airport  is  serving  a  mixed  type  of  aircraft  traffic,  all 
landings  could  possibly  be  converted  into  an  equivalent  aircraft.  One 
of  the  major  assumptions  of  the  analysis  procedure  was  that  only  the 
permanent  deflections  made  directly  under  the  wheel  were  considered. 

With  further  refinement  of  the  procedure,  the  permanent  deflections  away 
from  the  wheel  could  be  accounted  for. 

Procedure  No,  2 

The  second  solution  procedure  involved  the  initial  strain  concept 
commonly  used  in  the  finite  element  method.  A  brief  description  of  the 
initial  strain  concept  is  presented  below.  The  procedure  for  computing 
residual  strains  using  initial  strains  is  then  presented. 

According  to  Desai  and  Abel  (1972),  initial  strains  in  structural 
systems  may  occur  due  to  a  number  of  various  influences  as  change  in 
temperatures,  creep,  shrinkage,  crystal  growth  and  local  insitu  temper¬ 
ature  changes.  If  (£q)  is  a  vector  containing  the  initial  strains  and 
{e}  is  a  vector  containing  the  total  strains,  the  vector  {ee>  called 
the  effective  elastic  strain  vector  can  be  defined  by 

{t  }  -  {e}  -  {e  }  (19) 

e  o 

For  elastic  stress  strain  relationships 

(o>  -  [cKe  }  -  (c]{{e}  -  (e  }} 
e  o 


(20) 
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where  la}  is  the  stress  vector  and  [ c J  contains  the  material  elastic 
constants  that  relate  stress  and  strain.  The  strain  energy  function  can 
Chen  be  written  as 


U  =  j///v{e}TCc){e}dv  -  ///v{e}T(c){e0}dv 
+  I///v{eo}T(c){£o}dv 


(21) 


The  net  effect  on  equilibrium  from  the  initial  strains  is  an 
additional  load  term  {QQ}which  represent  the  contributions  of  the 
initial 
strains. 


[K] {q}  -  {Q}  +  {Qq}  (22) 

where  (Q  }  ■  fff  [B]T[c]{e  }dv 

o  v  o 

[k]  is  the  element  stiffness  matrix,  { q)  is  the  nodal  displacement 
vector,  {Q}  is  a  load  vector,  and  [B]  is  a  matrix  that  relates  strains 
and  displacements.  By  using  the  additional  load  vector  {Qq},  the 
effects  of  the  initial  strains  can  be  included  in  the  finite  element 
formulation.  The  initial  strain  procedure  above  was  termed  using 
elastic  material  laws. 

For  materials  that  behave  in  an  elastic  plastic  fashion,  the  method 

described  below  uses  the  same  stiffness  throughout  and  thus  shortens 

computational  time.  See  Figure  16.  Under  the  effects  of  load  {Q}  and 

elastic  stiffness  [k],  the  solved  displacement  will  correspond  to  point 

A.  The  correct  equilibrium  state  would  correspond  to  point  B  on  the 

non-linear  curve.  The  difference  between  the  correct  displacement 

(q,  }  and  the  computed  displacement  {q  }  is  {q  }.  The  "initial  strain", 

a  o 


i 


ess 


Displacements 


{V  {qb} 

Figure  16.  Initial  Strains  (after  Desca  and  Abel  (1972)) 


Figure  17.  Initial  Strain  Procedure  for  Repeated  Loadings  of 
Non-Linear  Materials 
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{eq}  corresponding  to  {qQ}  is  {eq}  ■  [B]  (qQ}.  The  correction  load  is 
found  by  using 

{Q  }  -  ///  [B][c]{e  }dv  (23) 

O  V  o 

The  load  is  applied  to  correct  the  displacements  from  A  to  B. 

Figure  17  represents  a  hypothetical  stress  strain  for  a  non-linear 

material  subjected  to  repeated  loading.  The  initial  strain  procedure 

described  above  could  be  used  to  adjust  the  initial  loading  to  point  B. 

The  strain  state  after  repeated  loading  would  correspond  to  Ec  and  the 

"initial  strains"  would  be  e  =  e  -  e_.  The  corresponding  correction 

o  c  B 

load  would  be 

{Q  >  -  ///  [B] [c] { {e}  -  { e_ }dv  (24) 

O  V  CD 

To  find  the  permanent  set,  S,  the  load  would  be  removed  using  the 
unloading  modulus. 

An  initial  attempt  was  made  to  employ  the  initial  strain  procedure 
for  computing  residual  displacements.  In  applying  the  procedure, 
tensile  strains  were  generated  in  certain  elements.  If  the  corres¬ 
ponding  tensile  stresses  exceed  the  tensile  strength  of  the  material, 
these  stresses  could  be  removed  by  reapplying  the  initial  strain  pro¬ 
cedure  or  by  eliminating  the  element  from  the  global  stiffness  matrix 
and  redistributing  the  load. 

The  development  of  a  complete  software  to  correctly  perform  this 
analysis  was  deemed  beyond  the  scope  of  this  feasibility  study. 

Conclusions 

In  this  chapter,  two  methodologies  for  computing  residual  displace¬ 
ments  in  surfaced  and  unsurfaced  pavement  systems  were  presented.  The 
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first  procedure  was  by  far  the  simplest  and  easiest  to  implement.  The 
procedure  requires  the  use  of  a  finito  element  computer  program  and 
requires  five  steps: 

(1)  Computing  element  strains  due  to  the  appropriate  wheel  loading  by 
finite  element  program. 

(2)  Computing  element  strains  due  to  unloading  conditions  by  finite 
element  program. 

(3)  Computing  the  permanete  set  strains  after  one  loading,  e^,  for  each 
element  by  subtracting  (2)  from  (1). 

(4)  Computing  the  permanete  set  strains  for  N  loadings  based  on  equa¬ 
tions  developed  from  laboratory  triaxial  testing. 

(5)  Computing  surface  deflections  by  summing  the  change  of  height  of 
each  element. 

No.  of  Elements 

S  -  £  (e  *h.) 

i-1 

Example  calculations  were  worked  using  the  first  procedure.  The 
first  two  cases  consisted  of  computing  residual  displacement  for  air¬ 
craft  tires  on  clay  and  an  expedient  crater  repair.  The  sample  calcu¬ 
lations  compare  favorably  with  published  data  on  the  actual  behavior. 

The  third  example  considers  a  possible  application  to  a  concrete  pave¬ 
ment  and  takes  into  account  the  fact  that  aircraft  traffic  is  distri¬ 
buted  across  the  runway's  cross  section. 

The  second  procedure  involves  the  application  of  the  initial  strain 
concept  in  finite  element  method  of  analysis.  The  initial  strain 
concept  would  be  used  to  formulate  corrective  nodal  loads  that  would  be 
applied  to  produce  permanete  displacements.  The  procedure  was  not 


easily  employed  In  that  initial  attempts  to  implement  it  produced 
tensile  strains  in  the  pavement  materials.  These  tensile  strains  would 
not  actually  exist  and  would  have  to  be  eliminated  from  the  finite 
elements  before  the  corrective  forces  could  be  applied. 
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CHAPTER  4 

Summary  and  Recommendations 


Summary 

Based  on  the  literature  survey  and  the  sample  calculations  in 
Chapter  3,  the  computations  of  residual  displacements  are  very  feasible. 
Two  different  methodologies  were  proposed.  The  method  that  was  the 
easiest  to  implement  requires  a  finite  element  computer  code  and  the 
results  of  laboratory  repeated  load  tests  on  pavement  materials. 

Finite  element  programs  were  used  to  obtain  the  residual  dis¬ 
placements  after  one  load  application.  The  use  of  non-linear  material 
behavioral  laws  or  were  recommended.  Relationships  between  the  first 
cycle  set  strain  and  the  set  strain  after  N  cycles  were  obtained  from 
the  triaxial  test  results.  These  relationship  produced  the  final  set 
strain  for  each  element,  e^.  The  surface  deflections  were  computed  by 
summing  the  settlements  of  each  element  or 


S  =  T.  eih1  (25) 

where  S  was  the  surface  settlement  and  h^  was  the  height  of  each 
element.  Sample  calculations  were  presented  to  show  the  effectiveness 
of  the  procedure. 


Recommendations 

The  recommendations  for  further  research  and  development  are  listed 
below: 

(1)  Develop  computer  software  to  implement  the  proposed  procedures. 

(2)  Conduct  extensive  laboratory  repeated  load  tests  on  a  wide  variety 
of  soils  and  pavement  materials  and  derive  general  repeated  load 


characteristics  that  would  be  related  to  fundamental  properties  as 
density,  void  ratio,  unconfined  compressive  strength,  moisture 
content  or  cone  penetration  resistance. 


(3)  Perform  repeated  load  tests  on  actual  pavements  or  crater  repairs 
where  the  induced  rut  depth  are  measured  throughout  the  testing 
program. 

(4)  Develop  software  to  implement  the  initial  strain  procedure  for 
computing  residual  displacements. 

(5)  Compare  which  constitutive  laws,  e.g.  Hardin,  overburden  dependent, 
linear,  or  hyperbolic,  yield  the  best  estimate  for  the  first  cycle 
permanent  strains. 
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